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MODIFIED  IMAGE  THEORY  QUASI-STATIC  RANGE  SUBSURFACE- 
TO-SUBSURFACE  AND  SUBSURFACE- TO- AIR  PROPAGATION  EQUATIONS 


INTRODUCTION 


During  the  past  few  years,  there  has  been  considerable  interest  in 
the  determination  of  the  quasi-static  field  components  of  antennas  loca- 
ted above  or  buried  within  the  earth's  surface,  where  the  quasi -static 
range  is  defined  as  the  range  where  the  measurement  distance  is  much 
less  than  a free-space  wavelength.  Quasi-static  range  results  are  use- 
ful for  submarine  radio  communication  and  detection  purposes,  as  well 
as  for  the  buried-miner  problem.  They  are  also  helpful  to  geophysicists 
engaged  in  determining  the  electrical  properties  of  the  earth. 

The  fields  from  various  subsurface  sources  previously  have  been 

derived  for  many  cases^-^  when  the  measurement  distance  R (=Vp2  ♦ Z2) 
is  much  larger  than  the  earth-skin  depth,  6 (i.e.,  the  quasi-near  range 
case),  and  when  R <<  6 and  |z  + h|  <<  5 (i.e.,  the  dc  case). 

Some  work  has  been  done^»b"^  on  determining  the  fields  from  vari- 
ous subsurface  sources  for  the  quasi -static  range  where  R is  comparable 
to  6.  However,  the  resulting  field-strength  expressions  are  very  com- 
plicated since  they  involve  products  of  modified  Bessel  functions  of 
different  argument.  Computing  the  field-component  expressions  is  a 
lengthy  and  difficult  process,  but  some  numerical  results  have  been 
obtained.  One  method  of  obtaining  these  results  has  been  discussed  by 
Atzinger,  Pensa,  and  Pigott.11  Numerical  integration  techniques  also 
have  been  employed. *0 

In  order  to  provide  simple  engineering  expressions  for  the  general 
quasi-static  electromagnetic  fields  produced  by  subsurface  antennas, 
finitely  conducting  earth-image  theory  techniques  will  be  employed. 

These  techniques  have  already  proved  quite  useful  in  determining  the 
quasi-static  fields  of  antennas  located  near  the  earth's  surface  — for 
both  singlelayered  and  multilayered  earths. 12-16 

Physically,  the  essence  of  the  finitely  conducting  earth-image 
theory  technique  is  to  replace  the  finitely  conducting  earth  by  a per- 
fectly conducting  earth  located  at  the  (complex)  depth  d/2,  where  d = 

2/y  = 6 (1  - i).  Analytically,  this  corresponds  to  replacing  the  alge- 
braic reflection  coefficient,  (u  - X)/(u  + X),  in  the  exact  integral 
equations  by  exp(-Xd),  where  X is  the  variable  of  integration. 
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For  antennas  located  at  or  above  the  earth's  surface,  the  gen- 
eral image-theory  approximation  is  valid  throughout  the  quasi-static 
range. 12,16  However,  if  one  or  both  antennas  are  buried,  the  previously 
derived1’  surface-to-surface  and  surface-to-air  image-theory  results 
(multiplied  by  the  exponential  attenuation  with  depth  factor  exp[y  (z  + 
h) ] or  exp(yh))  generally  will  be  valid  only  for  R > 3 |z  + h|,  where  h 
and  z,  respectively,  are  the  depth  of  the  source  and  receiving  antennas. 

Therefore,  we  will  modify  further  the  exact  subsurface-to-subsur- 
face  and  subsurface-to-air  integral  equations  and  obtain  fairly  simple 
field-strength  equations  for  the  general  quasi-static  range — which 
reduce  to  previously  derived  results  when  R >>  6 and  when  R <<  6 and 
|z  + h|  <<  6 . 

For  the  purposes  of  this  report  we  will  consider  six  buried  sources 
the  vertical  electric  dipole  (VED) ; vertical  magnetic  dipole  (VMD) ; 
horizontal  electric  dipole  (HED) ; horizontal  magnetic  dipole  (HMD); 
finite-length  horizontal  electric  antenna;  and  long  horizontal  line 
source.  All  these  sources  are  located  at  depth  h(h  <_  O')  with 
respect  to  a cylindrical  coordinate  system  (o,  4>,  z)  and  are  assumed  to 
carry  a constant  current  I.  The  VED  and  HED  (of  infinitesimal  length, 

£)  are  oriented  in  the  z and  x directions,  respectively,  and  the  axes 
of  the  VMD  and  HMD  (of  infinitesimal  area,  A)  are  oriented  in  the  z and 
y directions,  respectively.  The  earth  occupies  the  lower  half-space 
(z  < 0),  and  the  air  occupies  the  upper  half-space  (z  > 0).  Displace- 
ment currents  are  neglected  in  the  ground  and  (for  most  cases)  in  the 
air.  The  magnetic  permeability  of  the  earth  is  assumed  to  equal  pD, 
the  permeability  of  free  space.  Meter-kilogram-second  (MKS)  units  are 
employed  and  a suppressed-time  factor  of  exp(io)t)  is  assumed. 


SUBSURFACE-TO-AIR  PROPAGATION 


VERTICAL  MAGNETIC  DIPOLE  (VMD) 

When  h <_  0 and  z >_  0,  the  quasi-static  range  integral  expression 
for  the  VMD  Hertz  vector  is 


n 


z 


% 


2A 

u + A 


e'Az+uhJ0(Xp)dX  • 


The  fields  in  the  air  are  given 

dnz  ^nz 

L4>  " luvlo  aP  * Hp  = aPaz 


by 


(1) 


C2) 
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where 


yQ  = iwO0e0)  'v  0 (air), 

u 

y ^ (iwy0o)  (earth), 
u = (X2  + yw)  , and 

JQ  (Xp)  is  the  Bessel  function  of  the  first  kind,  order  zero,  and 


argument  Xp. 

Utilizing  the  Wait  and  Spies 


image 


theory  approximation, 


14 


u - X 

^ 

u + X 


-Xd 

e 


where  d = 6 (1  - i),  yields 


2X 

u + X 


% 


-Xd 

e 


Thus,  equation  (1)  becomes 


n ^ 


(1 


-Xd.  -Xz+uh.  , . . 
e )e  J0 ( Xp) dX 


(3) 


(4) 


(5) 


This  expression  for  flz  cannot  be  evaluated  analytically  (when  h 4 0) 
throughout  the  quasi-static  range.  However,  if  exp(uh)  can  be  replaced 
by  exp(f(Xh)]  times  exp(f(yh)j,  this  integral  can  be  evaluated  readily. 

Now  we  know  that  when  R >>  6 and  R >>  |h|,  exp(uh)  % exp(yh). 
Furthermore,  if  R <<  6 and  |h|  <<  6,  exp(uh)  'v  exp(Xh).  An  approxima- 
tion that  is  good  to  within  4 percent  if  A > B is 

V a2"  + 82  'v  0.96A  + 0.4B  . (6) 

Therefore,  we  will  approximate  exp(uh)  by 

uh  . yah  Xbh 

e ^ e e , (" ) 

with 

a = 0 and  b = 1 for  R <<  6 and  |h|  <<  6, 
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a = 0.4  and  b = 0.96  for  R/6  less  than  approximately  1, 
a = 0.96  and  b = 0.4  for  R/6  between  approximately  1 and  10,  and 
a = 1.0  and  b = 0 for  R > | 3h | . 

Thus,  equation  (5)  reduces  to 

oo 

„ IAeYah  / ,,  -Ad.  -A(z-bh)  ,,  ... 
n2  -v  — ^ — I (1  - e )e  v JJ0(Ap)dA  . 

0 

c-  1? 

Since 


then 


where 


(Ap)dA  = — ■■■■  - 
\P2  + S2 


n 

z 


Ki  = p2  + (z  - bh)2  and  = P2  + (d  + z - bh)2 


(9) 


(10) 


But  this  expression  is  identical  to  the  previously  derived  quasi - 
static  range  air-to-air  propagation  VMD  Hertz  potential  given  in  a 
previous  report, 12  equation  (36),  except  that 


1.  z - h is  replaced  by  z - bh, 

2.  z + h is  replaced  by  z - bh,  and 

Y dll 

3.  The  equation  is  multiplied  by  e 

Therefore,  the  subsurface-to-surface  and  subsurface-to-air  field- 
component  expressions  produced  by  a buried  VMD  (or  any  other  buried 
source)  may  be  obtained  from  those  given  in  the  previous  report!2  sim- 
ply by 

1.  Replacing  (z  - h)  by  (z  - bh), 

2.  Replacing  (z  + h)  by  (z  - bh) , and 


4 


■ 
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Y ah  - 

3.  Multiplying  each  equation  by  e (note  that  h ^ 0 ). 

Thus,  for  the  VMD  subsurface-to-air  propagation  case  (and  the  sub- 
surface-to-surface  case  — i.e.,  when  z = 0+) , equations  (37)  to  (39)  of 
the  previous  report^  become 


E.  - - ^ 


yah 


1 1 

JTT-kJ 


H 1A  e'fah(3o) 
p 4tt 


(z  - bh)  (d  + z - bh) 
K5  K5 


(ID 

(12) 


and 


H -v  - 
z 


IAeYah 

1 [i 

4tt 

K3 

- bh)2 


K3 

2 


3(d  + z - bh) 
1 - k2 

2 


1 


(13) 


VERTICAL  ELECTRIC  DIPOLE  (VED) 

For  the  VED  subsurface-to-air  propagation  case,  the  equations  of 
the  previous  report^  (equations  (29)  to  (31))  must  also  be  multiplied 
by  Yo2/^2-  Thus , 


2iro  \Ki  0 / 


n ^ 

z 


E -v 
P 


, . yah 
Hoe 

6 (z  - bh)  ya 

(d  + z - bh)  lz_ 

4ttd 

K5  p7 

k2 

ir1]) 


_ . Keyah  1 

z v 2tia  K3 


! 3(z  - 

K? 


- bh)2 

J * 


(14) 

(15) 

(16) 


and 


H.  -v 


He 


yah 


2ir 


(17) 


HORIZONTAL  ELECTRIC  DIPOLE  (HED) 

The  HED  equations  can  be  derived  from  equations  (4) , (7) , (8) , and 
(10)  to  (15)  of  the  previous  report. 12  They  are 


5 


^ 


I ft  COS  4*eyah  I 6p  (z  - bh 


I6P(z  - bh)  _4_TU  - bh)  (d  + z - bh)]| 

Kf  pd2  [ *1  K2  JJ 

H Ig  sin  [ (d  + z - bh)  (z  - bh) 

p 411  I K23  ' K3 

1_  (d  z - bh)  (z  - bh)H 

P2  [ *2  “ Kl  JJ  ’ 


.(23) 


z - bh)  fz  - bh) 


H % - M cos  »eYah  (<*  + z - bh)  _ (z  - bh)  1 
^ 4ttp^  K2  Kj  | 


H „ I ftp  sin  4*eyah  1 1 1 


z 4"  (K?  k||  • 

HORIZONTAL  MAGNETIC  DIPOLE  (HMD) 

The  appropriate  equations  for  HMD  subsurface-to-air  propagation  can 
be  obtained  from  equations  (48)  to  (56)  of  the  previous  report,^  and 
src 
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If  Kj  <<  6 and  |ah|  <<  5,  the  modified-image  theory  subsurface-to- 
air  (and  subsurface-to-surface  — i.e.,  z = 0+)  VED,  VMD,  HED,  and  HMD 
equations  (equations  (10)  to  (35)  of  this  report)  reduce  to  the  dc  equa- 
tions (for  example,  see  page  3-4  of  Kraichman2).  Furthermore,  if  Kj  » 

6 and  R » |h|,  the  equations  reduce  to  the  quasi-near  range  results 
given  by  Bannister4  (see  also  tables  3.13  and  3.14  of  Kraichman^). 
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VMD,  VED,  HED.  AND  HMD  AIR-TO-SUBSURFACE  PROPAGATION 


Now  that  simplified  expressions  for  VMD,  VED,  HED,  and  HMD  subsur- 
face-to-air  (and  subsurface-to-surface)  propagation  have  been  derived, 
simplified  expressions  for  VMD,  VED,  HED,  and  HMD  air-to-subsurface  (and 
surface-to-subsurface)  propagation  can  be  derived  easily  by  utilizing 
the  reciprocity  theorem.  This  theorem  (applicable  to  dipoles  in  the 
presence  of  any  linear  media)  states  that  voltage  V2  induced  in  antenna 
2 by  current  Ij  of  antenna  1 is  the  same  as  voltage  Vj  induced  in 
antenna  1 by  an  identical  current  I2  flowing  in  antenna  2.  Further 
details  are  given  in  several  references.1®-22 

Application  of  this  theorem  results  in 

_HM,  , . . „VE  A r. 

E (z,h)  = iujp0H  cos  <(>  - (h,z)  , (36) 

z 0 i. 

where  (a, 6)  = (height  or  depth  of  transmitting  dipole,  height  or  depth 
of  receiving  dipole); 


„HM,  VM  . x A™  ,.  . 

H,  (z,h)  = - H sin  $ -yrr  (h,z)  , 

z p a. 


(37) 


HE 

HF  VE  j? 

E7  (z,h)  = - Ep  cos  4 -yg-  (h,z) 


(38) 


and 


..HE , ..  -1  _VM  . A l ,,  . 

H (z.h)  = sin  $ - (h,  z)  , 

z iwuo  41  A 


HHE(z,h)  =_L_E»M^(h,z)  , 

p uuuq  4>  A 


n"E(z,h) 


-1  phm  1 n,  o 

T (h»z)  » 


iwuo  P A 


(39) 

(40) 

(41) 


_VE  „VM  hHM  „HM  _HE  HE  _ VE  VM  HM  HM  HE  HE 

Ez  * Hz  * ”p  * % * Ep  ’ E*  (2’h)  " Ez  ' Hz  ' Hp  * H4  ' Ep  ' (h,z)- 

(42) 
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BURIED  FINITE-LENGTH  HORIZONTAL  ELECTRIC  ANTENNA 


The  equations  for  subsurface-to-air  and  subsurface-to-surface  propa- 
gation, utilizing  a buried  finite- length  horizontal  electric  antenna 
(which  is  oriented  in  the  x direction)  as  a source,  may  be  obtained  from 
equations  (73)  to  (81)  of  the  previous  report,^  and  are 


T \ah 
„ If 

n 'v  

x 4thu)Eo 


£n 


Ml 


(x  + L/2) 


K12  - U - L/2) 


- tn 


'21 


(x  + L/2) 


K22  - (x  - L/2) 


IeYah 

n -v  - - — ; 

z 4ttiu£o 


Jin 


:n  - (d  + z - bh) 


Ki2  - (d  + z - bh) 


- Hn 


K21  - (z  - bh) 


K22  - (z  - bh) 


II 


(43) 


(44) 


and 


where 


V- It  'v 


K?2 


IeY3h 

1 

1 

2iro 

K2i 

k22 

(d  ♦ 

z - 

bh)2. 

(d  ♦ 

z - 

bh)2. 

(z  - 

bh) 

2 , and 

(z  - 

bh) 

2 ^ 

(45) 


In  addition 
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iuu0Ie 


yah  | 


4 ti 

d2 

fcx  + L/2) 

2 

rH 

ro  CM 

1 

Kn  - (x  + L/2)' 


K12  - (x  - L/2) 


K3 

22 


■]!■ 
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K21  - (x  + L/2) 


K22  - (x  - L/2) 


(46) 
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2 no 

K3 
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22 
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E — iJy2  * (d  + z - bh^ 

x 2tt  > y2  + (z  . bh)2 


(d+z-bh) 


(z  - bh) 


y 2 it  I [y2  + (d  + z - bh)2]  [y2  + (z  - bh)-  | 


, (54) 


_ Iye'an 1 1 _ _ 

Hz  % 2 tt  [y2  + (d  + z - bh)2]  [ y 2 + (z  - bh)2] 


. (55) 


SUBSURFACE-TO-SUBSURFACE  PROPAGATION 


HORIZONTAL  ELECTRIC  DIPOLE  (HED) 

When  h ^ 0 and  z <_  0~ , the  quasi-static  range  integral  expressions 
for  the  HED  Hertz  vectors  are 


Hx  ^ c~yR°  e~>R 

Ci  ' R0  " Ri 


r-  * / (z2H)'uU*h,J»<^>d> 


— 'v  - cos 
L1 


• • /(*)••' 


,J1(Xp)d>  , 


where 


cos  <j>|^-r(l  + yRo)e  + >Rl)e 


♦ 2 J Xeu(z+h)J1(Xp)dx|  , 


= 1 1/  (4ito ) , 
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R2  = p2  + (z  - h)2»  and 

R2  = P2  + (z  + h)2  . 

The  fields  in  the  earth  are  given  by 

£ = - y2n  + grad  div  n , 

159) 

H = o curl  n • 

These  expressions  (equations  (56)  to  (58)  for  II  , II  , and  7*11  (and 
the  resulting  field  components)  can  be  evaluated  analytically  throughout 
the  quasi-static  range.  In  fact,  they  have  been  - see  Bannister  and 
Mart. ® However,  the  resulting  field-strength  expressions  are  very  com- 
plicated because  they  involve  products  of  modified  Bessel  functions  of 
different  argument. 

From  equation  (6),  we  note  that 

u(z+h)  ya(z+h)  Ab(z+h)  (60) 

e 've  e . 


whe  re 


a = 0 and  b = 1 for  Rj/6  <<  1, 

a = 0.4  and  b = 0.96  for  Rj/6  less  than  approximately  1, 
a = 0.96  and  b = 0.4  for  Rj/6  between  approximately  1 and  10,  and 
a = 1 and  b = 0 for  p > 3 1 z + h|  . 


If  we  substitute  equations  (3)  and  (60)  into  equations  (56)  to 
(58),  we  can  easily  evaluate  the  MED  Hertz  potentials.  The  resulting 
expressions  are 


«x 

Ci 


e'YR° 

Ro 


e-rRi 


eya(z+h) 


1_ 

K, 


1_ 


(61) 


a ya(z+h) 
cos  4>e 

d - b(z  + h) 

, b(z  ♦ h) 

0 

Ku 

(621 


and 
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7«n 

Ci 


^ - p cos 


e-*K° e^Rl 2c 


-3—  (1  + yR0) 3-  (1  + yKi)  * 

Rn  R? 


ya(z+h) 

71 


, (63) 


'I  ”3 

where  K3  = p2  + [ b ( z + h]2  and  K2  = p2  + [d  - b(z  + h)]2. 

The  field-strength  expressions  may  be  determined  easily  by  employing 
equations  (59)  and  (61 j to  (03) . Thus, 


E -v- 
P 


I 8.  COS  $ 

e'  yR0 

4tto 

k 

)\  - y2U  - 


E . 'v 


£lY^i 

\(— 

R? 

Lv*i 

2eya(z+h) 

+ 1 

K3 

18  sin  $ 

e‘YR° 

4iro 

Ro 

(1  + yRq I 


(1  + yRjjl  - y2(z  + h) 2 


h)' 


)■ 


(64) 


3b3(z  + h)2 

1 + b - - — tj — — + yab(z  + h) 


K1 


,-tRi 


2eya(z+h) 

T1 


E ^ 


1 + 


24  (-£) 


(65) 


1£  cos  41 

P(z  - h) 

4tto 

R5 

K0 

♦ PU  * h)  (3  ♦ 3yR]  ♦ y^e'^ 


R1 


(66) 


p 4tt 


. ILZJH  (,  ♦ yR0)e-YR° 

RC 


(67) 


+ e 


y a( z+h ) 


d - b(z  + h)  + b ( z ♦ h)  + d - b(z  * h) 


K^P' 


k3p4 
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H x - 


It  CG  $ 

(z  - h) 

4 IT 

Rd 

R!3 


( )e 


-YRi 


ya(z+h) 


d - b (z 


and 


,2  ; 

[ 

I Up  sin  4> 

e”YR0 

4tt 

d3 

Ro 

♦ IQ  b(z  h)  1 

K3  Jl  * 


(68) 


e^Rl 

(1  + yR0) — (1  ♦ yRj) 


(69) 


+ e 


ya(z+h) 


1_  _ J_ 

Kl 


If  either  z or  h = 0,  then  Rq  = R]  and  the  first  two  terms  of  the  E , 
E$,  and  Hz  components  (equations  (64),  (65),  and  (69))  cancel  each  other? 
Furthermore,  if  h = 0,  the  first  two  terms  of  the  Ez  and  components 
(equations  (66)  and  (68))  double,  whereas  if  z = 0,  they  cancel  each 
other. 


VERTICAL  MAGNETIC  DIPOLE  (VMD) 

By  following  the  same  procedure  outlined  in  the  derivation  of  the 
HED  subsurface-to-subsurface  field-component  expressions,  the  VMD  sub- 
surface-to-subsurface  field-component  expressions  may  be  determined 
easily.  They  are 


ituwoIA  (e'YR° 


,-YRl 


E* * - — - ^ hr- (1  * YRo)  - tt  (1  + YRi) 


+ ya(z+h) 


..  IA  , . 
Hp  * 47  (P) 


*1  ' ij])  ' 

— R5  h)  (3  ♦ 3YR0  ♦ y2R§)e'YR° 


. 3[d  - b(z  ♦ h) J ya(z+h) 


K5 

K2 


(70) 


(71) 


and 


14 


» 


TR  5647 


H * “ 

z 4tt 


,-yRi 


,-yRq 


1 - 


3(z  - h)2 


(1  + yRq)  + Y P 


2„2  ! 


I . 3(2  ♦ h)2 


"f 


(l  + yRi)  + y2p2 


(72) 


+ e 


Ya(z+h) 


ic1 


, 3b2 (z  + h)2 

1 5 — — 

K3 


1 

5 


3[d  - b(z  ♦ h)]' 

* 2 

Ku 


If  either  z or  h = 0,  then  R0  = Ri  and  the  first  two  terms  of  the 
and  Hz  components  cancel  each  other.  Moreover,  if  z = h,  the  first  term 
of  the  Hp  component  expression  is  equal  to  zero. 

VERTICAL  ELECTRIC  DIPOLE  (VED) 

The  VED  subsurface-to-subsurface  field-component  expressions  also 
may  be  determined  by  following  the  procedure  outlined  in  the  derivation 
of  the  HED  subsurface-to-subsurface  field-component  expressions.  Thus, 


E 


It 

p(z  - h) 

4tto 

Ko5 

— [3  ♦ 3yR0  ♦ Y2R^]e"YR° 


and 


E -v  - 


p(z  + i 

R1 

It 

e'YR° 

4 ma 

R 3 

Ko 

■"f 


(1  + yRq)  + YZP2 


2„2 


H*  ' T,  <»’( 


(1  + yRi)  + Y^P 


.-yRi 


e yR0 

rT—  (1  * yR0)  - (1  + yRi) 

K0  K1 


(73) 


(74) 


(75) 


If  either  z or  h = O',  then  R0  = Ri  and  Ez  -v  0.  Furthermore,  if 

z = 0,  the  E.  component  doubles,  whereas  if  h = 0 , E 'v  0 (to  the 
0(Y2))- 
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HORIZONTAL  MAGNETIC  DIPOLE  (HMD) 

Again,  by  following  the  procedure  outlined  for  other  radiators,  the 
HMD  subsurface-to-subsurface  propagation  equations  are 


iuu  IA  cos  it 

n n 0 

(z  - h) 

P 4ir 

r3 

K0 

gYa(z+h) 
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P 2 

k4 
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: + h) 

K3  J)  ’ 


(76) 
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(80) 
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IA  sin  4> 
4tt 


(P) 


"I— -■  (3  + 3yR0  ♦ Y2R^]e"YR° 


R5 

0 


+ 3 f d - bCz  + h~)  1 e 


ya(z+h) 


K! 


(81) 


If  either  z or  h = 0",  then  Rq  = Ri  and  the  first  two  terms  of  the 
HMD  Hp  and  components  cancel,  and  Ez  ^ 0.  Furthermore,  if  z = 0,  the 
first  two  terms  of  the  E0  equation  double,  while  if  h = 0 they  cancel 
each  other. 

If  K3  <<  6 and  |a(z  + h) | <<  6,  the  subsurface-to-subsurface  HED, 
HMD,  VED,  and  VMD  equations  reduce  to  the  dc  equations  (for  example, 
see  page  3-4  of  Kraichman-) . Furthermore,  if  p >>  5 and  p >>  |z  + hj, 
the  equations  reduce  to  the  quasi-near  range  results. 


DISCUSSION 


Note  that  we  have  defined  loosely  the  variables  a and  b.  That  is, 
we  let  a = 0.96  and  b = 0.4  for  R 3 / 5 between  approximately  1 and  10, 
and  a = 0.4  and  b = 0.96  for  R 3 / 5 less  than  approximately  1.  The  speci- 
fic crossover  point  for  each  field-strength  component  will  depend  not 
only  on  R^S  but  also  on  |z  * h | / 6 . 

Since  the  resulting  field-strength  formulas  can  be  calculated 
easily  on  a desk  top  calculator,  the  field  strengths  can  be  determined 
by  using  both  values  of  a and  b,  thus  numerically  determining  the  cross- 
over point.  Alternatively,  we  could  set  the  two  expressions  (involving 
different  values  of  a and  b)  equal  and  solve  for  the  crossover  point. 

For  example,  consider  the  integral 

00 

N'  = |y  = J eU(z+h)J0(Xp)dA  , (82) 

0 

where 

N = I0  (A)  K0  (B) , 

A 3 2 tRl  + (2  * h)  J • 

B = j [Ri  - (z  + h) ] , and 
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I0  and  K0  are  modified  Bessel  functions  of  different  argument. 
The  (complicated)  analytical  solution  to  N'  is 


N'  "(2)  IoCA)K1  CB)  + I j (A)  K0  (B)  - (A)  Kj  (B)  - I1(A)K0(B)] 

and  the  (simple  form)  modified-image  theory  solution  is 

ya(z+h) 


N'  -v 


(83) 


(84) 


The  crossover  point  may  be  found  by  setting  the  two  versions  of 
(84)  equal  to  each  other  and  solving  for  o/5.  That  is,  let 


6 j N ' 


,al (z+h) / 6 


ia2(z+h)/6 


[ (P / 5)2  + b“(z  + h)2/62]*5  " [(p/5)2  + b^(z  + h)2/S2]^ 


(85) 


Thus , 


where  aj  = 0.4,  bj 
(a2  - a!) . 


(Z  * h) 

b?eP  - b^ 

6 

! - ep 

(86) 


0.96,  a2  = 0.96,  b2  = 0.4,  and  p = a[(z  + h)/6]  * 


If  (z  + h)/5  = -1.0,  the  crossover  point  is  at  o/5  = 0.46. 

The  numerical  integration  result  for  the  normalized  amplitude  of 
the  function  6M'  is  plotted  in  figure  1 versus  p/6  for  the  case  where 
(z  h)/6  = -1.  The  normalization  factor  (0  dB)  is  0.65,  the  value  of 
the  function  at  p/6  = 0.1.  Also  presented  is  the  modified-image  theory 
result  (with  a = 0.4  and  b = 0.96  for  p/6  < 0.46,  and  a = 0.96  and  b = 
0.4  for  p/6  > 0.46).  Here  we  see  that  the  modified-image  theory  result 
is  within  1 dB  of  the  numerical-integration  result  over  the  complete 
range  of  p/6  (0.1  to  10.0). 

As  another  example,  consider  the  magnetic  fields  at  a height  of  one 
skin  depth  (z/6  = 1.0)  produced  by  a HED  buried  at  a depth  of  one  skin 
depth  (h/6  = -1.0).  The  numerical-integration  results  for  the  normalized 
amplitude  of  each  component  (H'),  where 
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4H62H 
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f sin 

♦1 

( cos 

<P  1 

(87) 


are  presented  in  figures  2 and  3.  The  normalization  factor  (0  dB)  is 
the  numerical-integration  value  of  each  component  at  p/6  = 0.1. 


The  modified-image  theory  results  (from  equations  (24)  to  (26))  also 
are  presented  in  figures  2 and  3.  From  figure  2,  we  see  that  the  cross- 
over point  for  this  particular  example  is  p/6  = 1.5  (R/6  = 1.8).  That 
is,  if  0.1  < p/6  < 1.5,  then  a = 0.4  and  b = 0.96,  whereas  if  1.5  < p/6 
< 10,  then  a = 0.96  and  b = 0.4.  A comparison  (figure  3)  of  the  compos- 
ite modified-image  theory  results  with  the  numerical-integration  results 
shows  that  the  modified-image  theory  calculations  are  within  approxi- 
mately 1 dB  of  the  numerical-integration  results  over  the  complete  range 
of  p/6  (0.1  to  10).  In  fact,  the  modified-image  theory  even  predicts 
the  H0  component  amplitude  dip  at  the  right  place  (p/6  ■v  1.5). 


CONCLUSION 


Simple  expressions  for  HED,  HMD,  VED,  and  VMD  quasi -static  range 
subsurface-to-subsurface  and  subsurface-to-air  propagation  have  been 
derived  by  employing  finitely  conducting  earth-image  theory  techniques. 

It  has  also  been  shown  that  the  resulting  field-strength  expressions 
for  subsurface-to-air  propagation  are  the  result  of  simple  modifications 
to  one  of  the  author's  previously  derived^  quasi-static  range  air-to- 
air  propagation  equations. 

These  results  will  be  useful  for  submarine  radio-communication  and 
detection  purposes,  as  well  as  for  the  buried-miner  problem.  They  also 
may  be  helpful  to  geophysicists  engaged  in  determining  the  electrical 
properties  of  the  earth. 

Although  displacement  currents  have  been  ignored  in  the  analysis, 
they  can  be  included  simply  by  replacing  a by  o + iu>c  in  the  field- 
strength  equations  (as  long  as  |y  ( >>  jy^|). 

23 

In  a future  report  we  will  present  detailed  numerical  calculations 
of  the  simplified  field-strength  expressions  derived  in  this  report. 

They  will  be  compared  with  the  more  exact  numerical-integration  results 
whenever  possible. 
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Figure  3.  Comparison  of  Composite  Modified- Image  Theory 
and  Numerical- Integration  Results  for  the  Magnetic  Fields 
in  Air  Produced  by  a Buried  HED  (h/6  = -1.0,  z/6  = 1.0) 
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